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I. INTRODUCTION 
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Pseudospin symmetry is a relativistic symmetry of the Dirac Hamiltonian with scalar and vector 
mean fields equal and opposite in sign. This symmetry imposes constraints on the Dirac eigen- 
functions. We examine extensively the Dirac eigenfunctions of realistic relativistic mean field calcu- 
lations of deformed nuclei to determine if these eigenfunctions satisfy these pseudospin symmetry 
constraints. 
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Pseudospin doublets were introduced more than thirty years ago into nuclear physics to accommodate an observed 
near degeneracy of certain normal parity shell model orbitalswith non-relativistic quantum numbers (n r , I, j = £+1/2) 
and (n r — 1,£ + 2, j = £ + 3/2) where n r , £, and j are the single-nucleon radial, orbital, and total angular momentum 
quantum numbers, respectively 0,0- The doublet structure is expressed in terms of a "pseudo" orbital angular 
momentum, which is an average of the orbital angular momentum of the two orbits in the doublet, £ = £ + 1, coupled 
to a "pseudo" spin, s = 1/2 with j — £ ± s. For example, the shell model orbitals (n r Si/ 2 , (n r — l)d 3 / 2 ) will have 

£ = 1, (n r p 3 /2, (n r — l)/5/2) will have I = 2, for the two states in the doublet. Then the single-particle energy is 
approximately independent of the orientation of the pseudospin leading to an approximate pseudospin symmetry. 
These doublets persist for deformed nuclei as well Q. The axially-symmetric deformed single-particle orbits with 
non-relativistic asymptotic quantum numbers [N, n^, A]f2 = A + 1/2 and [N,n,3, A' = A + 2]S1' = A + 3/2 are quasi- 
& ' degenerate. Here N is the total harmonic oscillator quantum number, n 3 is the number of quanta for oscillations 
along the symmetry axis, taken to be in the z-direction, A and are respectively the components of the orbital and 
' q ■ total angular momentum projected along the symmetry axis Q|. In this case, the doublet structure is expressed in 
terms of a "pseudo" orbital angular momentum projection, A = A + 1, which is added to a "pseudo" spin projection, 
jl = ±1/2 to yield the above mentioned doublet of states with f2 = A — 1/2 and Q' = A + 1/2. This approximate 
pseudospin "symmetry" has been used to explain features of deformed nuclei, including superdeformation [f| and 
j ■ identical bands [ESIEI3 as well. 

Although there have been attempts to understand the origin of this "symm etry " [Hi ITl| . only recently has it 
been shown to arise from a relativistic symmetry of the Dirac Hamiltonia n Il2l Il3| which we review in Section II. 
This relativistic symmetry implies conditions on the Dirac eigenfunctions 1141 which we discuss in Sections II and 
III. These relationships have been studied extensively [Pj, [l5|, [la, Il7l Il8l Il9| for spherical nuclei. For deformed 
nuclei, the relationships have been studied only in a limited way and primarily for the lower components of the Dirac 
eigenfunctions [2(], l2ll |22^| . In this paper we shall test thoroughly these relationships between the upper and lower 
components of the two states in the doublet for realistic deformed relativistic eigenfunctions [2^, 0] • 
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II. THE DIRAC HAMILTONIAN AND PSEUDOSPIN SYMMETRY 



The Dirac Hamiltonian, H, with an external scalar, Vs(r), and vector, Vv(r), potentials is given by: 

H = &-p + p(M + V s (r)) + V v (r) , (1) 

where a, (3 are the usual Dirac matrices, M is the nucleon mass, and we have set Ti = c = 1. The Dirac Hamiltonian 
is invariant under a SU(2) algebra for two limits: Vs(r) = Vv(r) + C s and Vs(r) — —Vv{r) + C ps where C s ,C ps are 
constants |2^| . The former limit has application to the spectrum of mesons for which the spin-orbit splitting is small 
26] and for the spectrum of an antinucleon in the mean- field of nucleons [2ll |28| . The latter limit leads to pseudospin 
symmetry in nuclei 01 ■ This symmetry occurs independent of the shape of the nucleus: spherical, axial deformed, 
or triaxial. 



A. Pseudospin Symmetry Generators 

The generators for the pseudospin SU(2) algebra, Si (i = x,y,z), which commute with the Dirac Hamiltonian, 
[H ps , Si] =0, for the pseudospin symmetry limit Vg(r) — — Vv(r) + C ps , are given by 0] 

§i= (s, 0\(U p stU p 0\ 



„ s l J \ 

where Si = Ci/2 are the usual spin generators, a L the Pauli matrices, and U p = 17 p — is the momentum- helicity 

unitary operator introduced in Thus the operators Si generate an SU(2) invariant symmetry of H ps . Therefore, 
each eigenstate of the Dirac Hamiltonian has a partner with the same energy, 

where k are the other quantum numbers and // = ±^ is the eigenvalue of S z , 



The eigenstates in the doublet will be connected by the generators S± = S x ± iS, 



The fact that Dirac eigenfunctions belong to the spinor representation of the pseudospin SU(2), as given in Eqs. Q)- 
(JSJ), leads to the conditions on the corresponding Dirac amplitudes that are explored in this paper and developed in 
the next Subsection. 



B. Dirac Eigenfunctions and Pseudospin Symmetry 



An eigenstate <J>| S _ (f) of the Dirac Hamiltonian H ps , Eq. ©, is a four-dimensional vector, 



(r) 



(6) 



where g~ . (r) are the "upper Dirac components" and f~ . (r) are the "lower Dirac components" . The superscript 
+ (— ) indicates spin up (spin down). 
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The connections between the Dirac eigenstates of the doublet (jj, = ±i) resulting from Eqs. 0-101 lead to 
relationships between the Dirac amplitudes in Eq. 14], 

ft i(f) = /t"i(* ! ) = 0. (7a) 



f£Ar) = = , (71.) 



re, 2 re, 2 

ffjx(r) = -<^_i(r> = <?re(r> - ( 7c ) 



and to first order differential equations, 



d d \ _ ( d d \ , ,„ s 

fe"^J 5 ^ (r) = l^ + ^JC^' (8a) 
l^ = ± (£ Ti l)^- (8b) 

Thus pseudospin symmetry reduces the eight amplitudes for the two states in the doublet to four amplitudes. In the 
next Section we shall discuss these relations in Eqs. 0-® for axially deformed nuclei. 

III. PSEUDOSPIN SYMMETRY FOR AXIALLY DEFORMED NUCLEI 

If the potentials are axially symmetric, that is, independent of the azimuthal angle (j>, Vs,v{r) = Vs,v{p, z), p = 
\/ x 2 + y 2 , then the Dirac Hamiltonian has an additional U(l) symmetry in the pseudospin limit. The conserved U(l) 
generator is given by 

where l z = U P £ Z U P and £ z = (r x p) z . In this case, the Dirac eigenstates of H ps are simultaneous eigenstates of L z , 
with eigenvalue A, and of the total angular momentum generator J z = S z + L z , with eigenvalues Q = A + jl = A±^ 

Here fj denotes additional quantum numbers that may be needed to specify the states uniquely. 

The conventional method of labeling the eigenstates of axially deformed single-particle states in nuclei is to use 
the asymptotic quantum numbers (N, n.3, A, O), mentioned in the Introduction, that emerge in the limit of a non- 
relativistic axially-symmetric deformed harmonic oscillator with spin symmetry. For the relativistic axially-deformed 
harmonic oscillator with spin symmetry |29j| the eigenfunctions can also be labeled by these quantum numbers. 
However, only the spatial amplitudes of the upper components of the doublet will necessarily have the nodes suggested 
by these quantum numbers, whereas the spatial amplitudes of the lower components may have different nodal structure. 
For spherically symmetric potentials a general theorem relates the nodal structure of the upper and lower Dirac 
amplitudes, and has been used to explain the non-relativistic radial quantum numbers characterizing pseudospin 
doublets in spherical nuclei 30] . A corresponding theorem for axially-deformed potentials in the pseudospin and spin 
limits of the Dirac Hamiltonian appears to hold under certain conditions which the relativistic harmonic oscillator 
satisfies, but which do not generally apply for realistic axially-symmetric potentials (U- For the latter, only the 
quantum numbers A and fl in Eq. (|10|l are conserved in the pseudospin limit. The fact that the axial-symmetry of the 
potentials determines the (^-dependence of the Dirac wave functions, leads to the following form for the relativistic 
pseudospin doublet eigenstates [l4| 



-gl A (p,z)e^ 







n = A-i, (lla) 



4 



V 



n' 



A + i 



(lib) 



/ 



The two states in the doublet have the same pseudo-orbital angular momentum projection along the symmetry 
axis, A, but different total angular momentum projections il = A — | and O' = A + |. As seen from Eqs. (|lla|) and 

(|llb|l . the pseudospin projection, p — ±i, and A are respectively the ordinary spin projection and ordinary orbital 
angular momentum projection of the non- vanishing lower component /- ^(p, z). The corresponding dominant upper 

components <7- . _i (p, z) and gZ . i (p, z) have orbital angular momentum projections A = A — 1 and A' = A + 1 

?7,A, 2 77sA 2 

respectively, hence A' = A + 2. Accordingly, il = A + 1/2 and Of = A' — 1/2 = A+ 3/2. These assignments agree with 
the non-relativistic pseudospin quantum numbers discussed in the Introduction. The generic label fj in & p ~ s ^ . ^(r) 

replaces the harmonic oscillator labels N and rz.3, which are not conserved for realistic axially-deformed potentials in 
nuclei. 

In obtaining the expressions in Eq. l|llf) we have used the relations in Eq. Q , which for axially-deformed potentials 
read Q 



iO, z ) = ff,,A(P> Z ) 
'»7,a( 



(12a) 
(12b) 
(12c) 



and the differential relations 




d_± 

dz 9 vXt 



i{p,z) = ± 





A-l 


dp 


P 








•j) 


{ dp 





',A,± 



(13a) 
(13b) 



We shall now test to see if the pseudospin symmetry conditions in Eqs. H12fl - I|13|) are va lid for realistic relativistic 
mean field eigenfunctions in deformed nuclei. 



IV. COMPARISON WITH REALISTIC RELATIVISTIC EIGENFUNCTIONS 

The single-particle energies and wave functions for 168 Er are calculated by the relativistic Hartree theory with the 
parameter set NL3 in a Woods-Saxon basis [23L |24| . This method has been developed from relativistic theory in 
coordinate space [13, IH, and has the advantage that it easily generalizes to include both deformation and pairing 
correction self-consistently. The pairing correlation is treated with the BCS approximation. These calculations lead 
to a theoretical average binding energy B/A = 8.107 MeV, a quadrupole deformation (3 = 0.3497 and a root mean 
square radius R = 5.376 fm, which reproduce the data well. For these realistic eigenfunctions the harmonic oscillator 
quantum numbers are not conserved, but the orbitals are labeled by the quantum numbers of the main spherical basis 
state in the expansion of the dominant upper component in the Dirac eigenfunctions. 

In Figure 1, the calculated single-neutron energies, e — E — M, for the pseudospin doublets in 168 Er are presented. 
From left to right, the panels correspond to the pseudo orbital angular momentum projection A = 1, 2, 3 and 4, 
respectively. The energy splitting between members of pseudospin doublets decreases as the single-particle binding 
energy, -e, decreases. For pseudospin doublets with binding energy larger than 5 MeV, the spin-up (pseudospin down) 
state is higher than the spin-down (pseudospin up) one. On the other hand, for the bound doublets with binding 
energy less than 5 MeV, the opposite is observed. 

Four pairs of neutron pseudospin partners are chosen to illustrate the relations given above, (i) The states [402] | 
and [404] I (A = 3), which have a large energy splitting (about 2 MeV). The single-particle energies are respectively 
e [402]| = -12.083 and e [404] 7 = -14.160 MeV. (ii) The states [400]^ and [402] f (A = 1), which have a small energy 



5 



splitting (about 0.4 MeV). The single-particle energies are respectively £[4oo]§ = —10.2073 and £[404] J = —10.603 MeV. 
(iii) The states [501] | and [503] | (A = 2), which have a small energy splitting (less than 0.4 MeV). The single-particle 
energies are respectively emoija = —1-349 and £[ 50 3]5 = —0.9603 MeV. (iv) The states [510]^ and [512] | (A = 1), 
which have a tiny energy splitting (less than 0.01 MeV). The single-particle energies are respectively £[5io]± = —3.8436 
and £ [512] | = -3.8378 MeV. 

Plots for the above four pairs of neutron pseudospin partners are shown in Figures 2, 4, 6, 8, as a function of z for 
three segments: p = 1, 3, 5 fm, and in Figures 3, 5, 7, 9, as a function of p for three segments: z = 1, 3, 5 fm. In each 
segment the top row displays the relationships between lower component amplitudes given in Eqs. (|12afl - l|12b|) . and 
the relationship between upper component amplitudes given in Eq. (|12c|) . The bottom row displays the differential 
relationships between upper component amplitudes given in Eq. (|13fl . 

From these figures, we can draw a number of conclusions. First, while the amplitudes ft- 1 (p, z), f7r 1 (p, z) are 

not zero as predicted by Eq. (|12afl . they are much smaller than /J^ 1 (p, z), 1 (p,z). Furthermore, i{p,z) 

and f7~. _ L {p,z) have similar shapes as predicted by Eq. I|12bjl . In Figures 2-3 for the [402]5/2, [404]7/2, A = 3 

doublet there is some discrepancy in the shapes but the shapes become more equal as both the pseudo-orbital angular 
momentum projection A decreases (see Figures 4-5 for the [400]l/2, [402]3/2, A = 1 doublet) and the binding energy 
decreases (see Figures 6-7 for the [501]3/2, [503] 5/2, A = 2 doublet). 

The amplitude —gZ~. i{p,z) has the same shape as the amplitude gZ ~. 1 (p 7 z), in line with the prediction 

*7,A,— ^ ?7,A,^ 

of Eq. p2c(l . but they differ in magnitude. Again the discrepancy decreases as the pseudo-orbital angular 
momentum projection A decreases (compare Figures 2-3 for the [402]5/2, [404]7/2, A — 3 doublet with Figures 
8-9 for the [510]l/2, [512]3/2, A = 1 doublet) and the binding energy decreases (compare Figures 2-3 for the 
[402J5/2, [404J7/2, A = 3 doublet with Figure 6-7 for the [501]3/2, [503]5/2, A = 2 doublet or Figures 8-9 for the 
[510]l/2, [512]3/2, A — 1 doublet). These amplitudes are much smaller than the other upper amplitudes, gt- 1 (p, z). 

The differential relation in Eq. I)13a(l between the dominant upper components, gZ 1 {p,z) and gZ r x (p, z), is 

well obeyed in all cases. The differential relations in Eq. I|13b|l relate the dominant upper components, <?_ t 1 (p, z) 

to the small upper components g^ - 1 (p, z). The shapes of the left-hand-side and of the right- hand-side of Eq. (|13b[) 

are the same, but the corresponding amplitudes are quite different. Therefore, the differential relations in Eq. (|13b|) 
are less satisfied. These differences might partly originate from the differences in the magnitudes of the small upper 
components in Eq. I|12cfl . 

V. SUMMARY 

We have reviewed the conditions that pseudospin symmetry places on the Dirac eigenf unctions. We have shown that 
the conditions on the lower amplitudes, Eqs. H12a|) - (jl2b|) . are approximately satisfied for axially deformed nuclei. The 
differential relation between the dominant upper component amplitudes, Eq. (|13afl . is also approximately satisfied. 
However, both the relation between the amplitudes of the small upper components, Eq. I|12cfl . and the differential 
equations, Eq. i|13b[) . that relate the dominant upper components with the small upper components are not well 
satisfied. The pseudospin symmetry improves as the binding energy and pseudo-orbital angular momentum projection 
decrease, which is consistent with previous tests of pseudospin symmetry in spherical nuclei. 

This research was supported in part by the United States Department of Energy under contract W-7405-ENG-36, 
in part by a grant from the US-Israel Binational Science Foundation and in part by the Major State Basic Research 
Development Program Under Contract Number G2000077407 and the National Natural Science Foundation of China 
under Grant No. 10025522, 10221003, 10047001 and 19935030. 
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FIG. 3: Wave functions in (Fermi) 3 ^ 2 as a function of p and z = 1, 3, 5 fm for the neutron pseudospin doublet [402]5/2 and 
[404J7/2 (A = 3) in 168 Er. The content of the graphs in each segment as in Fig. 2. 
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FIG. 4: As in Fig. 2 but for the neutron pseudospin doublet [400]l/2 and [402]3/2 (A =1) in 168 Er. 
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FIG. 5: As in Fig. 3 but for the neutron pseudospin doublet [400]l/2 and [402]3/2 (A =1) in 168 Er. 
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FIG. 6: As in Fig. 2 but for the neutron pseudospin doublet [501]3/2 and [503]5/2 (A =2) in 168 Er. 
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FIG. 7: As in Fig. 3 but for the neutron pseudospin doublet [501)3/2 and [503]5/2 (A =2) in 168 Er. 
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FIG. 8: As in Fig. 2 but for the neutron pseudospin doublet [510]l/2 and [512]3/2 (A =1) in 168 Er. 
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FIG. 9: As in Fig. 3 but for the neutron pseudospin doublet [510]l/2 and [512]3/2 (A =1) in 168 Er. 



